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A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., 
incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The 
numerical method is based on a finite volume method that incorporates a pressure-staggered 
mesh and an incremental pressure equation for the conservation of mass. Comparison of three 
generally accepted time-advancing schemes, i.e.. Simplified Marker-and-Cell (SMAC), 
Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (IT A) scheme, 
are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over 
circular and square cylinders. Calculated results show that the ITA is the most stable 
numerically and yields the most accurate results. The SMAC is the most efficient 
computationally and is as stable as the ITA. It is shown that the PISO is the most weakly 
convergent and it exhibits an undesirable strong dependence on the time-step size. The 
degenerated numerical results obtained using the PISO is attributed to its second corrector step 
that cause the numerical results to deviate further from a divergence free velocity field. The 
accurate numerical results obtained using the ITA is attributed to its capability to resolve the 
nonlinearity of the Navier-Stokes equations. 

The present numerical method that incorporates the ITA is used to solve an unsteady 
transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a 
vitiated supersonic airstream. The turbulence fields in these flow cases are described using 
multiple-time-scale turbulence equations. 

For the unsteady transitional over an oscillating airfoil, the fluid flow is described using 
ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. 
It is shown that the numerical method successfully predicts the large dynamic stall vortex 
(DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating 
airfoil. The calculated streaklines are in very good comparison with the experimentally 
obtained smoke picture. The calculated turbulent viscosity contours show that the transition 
from laminar to turbulent state and the relaminarization occur widely in space as well as in 
time. The ensemble-averaged velocity profiles are also in good agreement with the measured 
data and the good comparison indicates that the numerical method as well as the multiple- 
time-scale turbulence equations successfully predict the unsteady transitional turbulence field. 

The chemical reactions for the hydrogen in the vitiated supersonic airstream are described 
using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used 
to describe the fine details (such as the instability) of chemically reacting flows while a 
reduced chemical kinetics can not be used confidently due to the uncertainty contained in the 
reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult 
to obtain a fully converged solution due to the coupling between the large number of flow, 
turbulence, and chemical equations. The numerical results obtained in the present study are in 
good agreement with the measured data. The good comparison is attributed to the numerical 
method that can yield strongly converged results for the reacting flow and to the use of the 
multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel 
and the oxidant. 



PAiu'c 



iNTEMMUY 5 IMt. 


1857 




iA 



< 

CO 

< 


1858 


* Resident Research Associate 



Introduction. 



CO 

<L> 

s- 


O 

© 

c 

o 

* 

+-> 

O 

oo 

£ 

E 

>> 

Tj 

cd 

© 
•4— > 

co 

c 

p 

<4-H 

o 

c 

o 

co 


Oh 

e 

6' 


a 'v 

g 

O 
© 

a 
© 

43 

O w 

^ & 
50 © 

c U 

• wi 

O T3 


<D 

> 

3 

© 


U 

< 

2 

00 


>v 

> 

T3 

c 
a 

H 1 H 


© 

I 


© 

ta 

43 


o 

00 


co 

© 

03 

Vh 

8 . 

o 


50 

C 

•3 

■4— > 

43 

Oh 

00 

•*-» 

• I— ( 

o 

43 

Oh 


a 

§ 

co 

© 


KH 00 


CO 

£ 

E 

^© 

*04 

I 

wi 


C/5 

a 


50 

J3 

3 

a 

C /5 

£ 

*§ 


£ 

o 


Oh 

g 


I 

T3 


© t_; 

T3 © 

o -© 

43 £ 

>> 3 

O >> 
H ° 

43 £ 
o a 
& a 

o ? 
cd cd 

Vh »-h 

© <D 


> 

O 


> 

O 


£ £ 

o o 

E E 



r> 



50 

C 

3 

33 

8 

O 


© 

> 

o 


£ 

o 


E 


13 

(3 




HH 

> 


1859 


Combustion of H2 in Vitiated Supersonic Airstream. 



Fluid flow inside Space Propulsion Systems includes all flow velocities, i.e., 
Incompressible Rows (M < e), 

Low Mach Number Rows, 
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Need to identify a best time-integration scheme among many available methods. 

Multiple-time-scale turbulence equations are used in the present unsteady 
transitional flow and the chemically reacting flow calculations. 
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1. Iteartive Time- Advancing Scheme (ITA) 
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2. Simplified Marker and Cell (SMAC) 
Predictor Step : Solve momentum equation. 
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<|) T 3d + i-ud = 



Solve eqs. (1-4) for each time-level.The velocity field strongly satisfy the 
conservation of mass. 
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3. Pressure-Implicit Splitting of Operators 
Predictor Step : Solve momentum equation. 
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Polar cavity flow starting from rest 
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Polar cavity flow starting from rest 
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Strouhal Numbers (St = 2fr/Uo) for Row over a Circular Cylinder 
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Ref. [2]: Braza, Chassing, Ha Minh, JFM 1986, 
Ref. [14]: Eaton, JFM 1987. 



Calculation of Unsteady Transitional Flow 
over Oscillating Airfoil* 
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Kim, Zaman, and Panda, To be presented in the AS ME Fluid 
Engineering Conference, Washington D.C., June 20-24, 1993. 
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Conservation of mass equation 
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Ensemble average 
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(i) The hot wire can not accurately measure the velocity components when the 
flow is misaligned more than approximately 30° from the hot wire axis. 

(ii) The interaction between the DSV and the TEV occurs in a relatively coarse 
mesh region and the numerical method yields somewhat deteriorated results. 











(g) 7.4° 

Turbulent viscosity contour ( = 10) 
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Chemical species concentration equation 
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Chemical reactions for the combustion of H2 in a vitiated supersonic airstream are 


3G Os 
- 00 
O ON 


<U cj 

T3 e3 


ft v 

o <§ 
c too 

q q 

cd X2 

0 0 

.2 >> 

1 3 

.S .a 

0 E 

a a 

fTt r*-i 


< o 


p c/3 

c3 -g 
<■> § 
«3 

.a SJ 

o E 
e - 

3 .2 

13 a 

o 3 


73 "" 4 

<D 73 
O <L> 

3 q 

"8 1 

u a 



r> 


£ 

3 

O 

M — 1 

cd 

q=! 


G 

O 

c3 

■4—4 

X) 

O 

Vh 

D 

X) 

O 

■ 4 — » 

S 

q 

■4—* 

q 

3 

<D 

0 

m _d 

W) 

^-H 

Ui 

C+H 

cd 

3 


■ 4 — • 

• 1 — ( 

<D 

45 

D 


44 

q 

cd 

<u 

g 

<L> 


£ 


■ 4 — > 



3 

X) 


too 

q 

q 

>^H 

C/3 


■ 

6 

3u 

q 

<u 

O 

45 

0 

0 

<D 

4— > 

(D 

45 

cd 

■ 4 —* 



<D 

O 

-4-* 


«) T 

C ^3 

O 

Vh ^ 
■D p 

CZ3 

C/3 

<u q 

^ .2 

2 § 

vs cr 

73 <u 


« q 
■o £ 
0 3 

45 -O 

s a 


5 JD 
q s 


£ J3 

H *-* 

C/3 O 

3 5/5 

O 73 

ti C 
ctf <D 


,S 9 ^ w 

^ ci 

'q S >> 
« q § O 
3 o 4: « 


O 1 2 

Cd C/3 C3 ^ 


o *0 <D tt 

a £? g g> 

a s -e c 

£ ° 3 o 
E— ' o q o 


1884 




CONCLUSIONS AND DISCUSSION 
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The ITA can best resolve the nonlinearity of the Navier-Stokes Equations. 
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conservation of mass, and 

(iii) the Multiple-time-scale turbulence equations that can resolve the 
inequilibrium turbulence field. 



